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Abstract 

This chapter provides an introduction to the fundamental physical 
ideas and models relevant to the phenomenon of magnetic hysteresis in 
nanoparticle assemblies. The concepts of single-domain particles and su- 
perparamagnetism are discussed. The mechanisms of magnetization by 
coherent rotation and the role of temperature in the gradual decay of 
magnetization are analyzed in the framework of simple analytical models. 
Modern numerical techniques (Monte Carlo simulations, Magnetization 
Dynamics) used to study dense nanoparticle assemblies are presented. An 
overview of the most common experimental techniques used to measure 
the magnetic hysteresis effect in nanoparticle assemblies are presented and 
the underlying principles are exposed. 
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1 Introduction 



Magnetic nanoparticles (MNPs) arc minute parts of magnetic materials with 
typical size well below 10 _7 m. They are present in different materials found 
in nature such as rocks, living organisms, ceramics and corrosion products, 
but they are also artificially made and used as the active component of fer- 
rofluids, permanent magnets, soft magnetic materials, biomedical materials and 
catalysts. Their diverse applications in geology, physics chemistry, biology and 
medicine renders the study of their properties of great importance both to sci- 
ence and technology. 

In geology, the nature and origin of magnetic phenomena related to the 
presence of magnetite nanoparticles in rocks is of great interest to the palaeo- 
magnctist who searches for the geomagnetic record of rocks. The presence of 
magnetite particles associated with the trigeminal nerve in pigeons offers a re- 
liable explanation to the Earth's magnetic field detection and the consequent 
navigation capability. In fine arts, magnetic analysis of ancient paintings facil- 
itates the reconstruction of the production techniques of ancient ceramics. In 
living organisms, the role of ferritin, a magnetic nanoparticle per se, is impor- 
tant among the iron storage proteins. MNPs are also used as contrast agents 
in Magnetic Resonance Imaging. Recent work has involved the development of 
bioconjugated MNPs, which facilitated specific targeting of these MRI probes 
to brain tumors. MNPs are also used as highly active catalysts which has long 
been demonstrated by the the use of finely divided metals in several reactions. 
Owing to their high surface-to- volume ratio MNPs of iron are more efficient at 
waste remediation than bulk iron. 

High density magnetic data storage media provide a major technological 
driving force for further exploration of MNPs. It is expected that if MNPs 
with diameter 5nm can be used as individually addressed magnetic bits, mag- 
netic data storage densities of 1Tbit /in 2 would be achieved, namely an order 
of magnitude higher than the present record (Moser 2002). MNPs have also 
been demonstrated to be functional elements in magneto-optical switches, sen- 
sors based on Giant Magneto-Resistance and magnetically controllable Single 
Electron Transistor devices. 

The most common preparation methods for MNPs produce assemblies with 
different structural and compositional characteristics that depend on the par- 
ticular method adopted. Granular films, fcrrofhiids and cluster-assembled films 
are characterized as assemblies with random order in MNP locations, while or- 
dered arrays are found in patterned media (known also as magnetic dots) and 
self-assembled films. The MNP preparation methods are divided to top-down 
and bottom-up. In top-down methods , the NPs are formed from a larger sys- 
tem by appropriate physical processing, such as thermal treatment, etching, 
etc. In bottom-up methods, the NPs are formed by an atomic nucleation pro- 
cess that takes place either in ultrahigh vacuum or in a liquid environment. 
The latter method relies on colloidal chemistry techniques and presently ap- 
pears to be the most promising method for production of nanoparticles with 
extremely narrow size distribution. Colloidal synthesis methods combined with 
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self-assembly methods produce MNP samples with both size uniformity and long 
range structural order. It is worth noticing that structural order in a MNP as- 
sembly is a decisive property for production of ultrahigh density storage media. 
Owing to their attractive features and their low cost, colloidal synthesis meth- 
ods and self-assembly attract presently intense research activity in the field of 
MNP preparation (Petit 1998, Murray 2001, Willard 2004, Farrcll 2005, Darling 
2005). 

The magnetic properties of MNPs and their assemblies provide a fascinating 
field for basic research, which is done on two different scales, the atomic and 
the mesoscopic. In the atomic scale, the properties of individual MNPs are 
examined and they are revealed in samples with low particle concentration. 
In the mesoscopic scale, dense samples are examined which exhibit collective 
magnetic behavior arising from interparticle interactions. The study of the 
magnetic properties can be naturally divided in the investigation of the ground 
state configuration (long range order, disorder, etc) and the excitations from it. 
Excitations can be either weak, as for example at low temperature and weak 
external magnetic field, or strong, as for example, close to a thermal phase 
transition or under a reversing magnetic field. 

For individual MNPs the ground state configuration can differ remarkably 
from the parent bulk material in various ways. For example, owing to energy 
balance reasons, the abundance of magnetic domains that form in a bulk magnet 
can be replaced by a single domain in a MNP, which then becomes magnetically 
saturated even in the absence of an external magnetic field (Neel 1949). The 
application of an external field forces the atomic magnetic moments of a single- 
domain MNP to rotate coherently (Stoner 1948). Also, for temperature above a 
threshold, the direction of particle's magnetization fluctuates at random, making 
the particle bahave as a molecule with a giant magnetic moment. The appli- 
cations of this effect, known as superparamagnctism (Bean 1959), are presently 
a lot, ranging from geology to medicine. Finally, we should remark that the 
above described simplified picture of a single-domain MNP becomes invalid if 
one considers the crucial effect of the MNP surface. Reduced crystal symmetry 
and chemical disorder close to the surface can produce variations between the 
surface and interior magnetic structure and modify the overall response of the 
MNP to an applied field (Kodama 1999). 

When MNPs form dense assemblies, interparticle interactions produce a col- 
lective behavior, by coupling the magnetic moments of individual MNPs. This 
fact renders in most cases even the determination of the ground state configura- 
tion an intricate physical problem. The collective behavior of dense (interacting) 
assemblies is reflected also on the modified magnetic response of the assembly, 
compared to isolated MNPs. The most complex behavior occurs in samples 
with random morphology and long-range magnetostatic interactions. Various 
experimental measurements have been proposed to reveal the nature of the in- 
terparticle interactions, and various measuring protocols probe different aspects 
of the collective behavior. On the other hand, analytical models have difficul- 
ties in predicting or explaining the magnetic behavior of these interacting MNP 
assemblies, and most of the curret research relies on numerical simulations. 
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In this chapter we provide an introduction to the fundamental ideas and 
concepts pertaining to the magnetic properties of MNP assemblies. Emphasis 
is given to the response of MNP assemblies to an applied magnetic field and the 
related issue of magnetization reversal. The chapter is organized as follows : 
In Section [5] we discuss the magnetic properties of individual (isolated) MNPs. 
Fist, the condition under which a single-domain MNP is formed is derived, and 
then the magnetic response under an applied field is examined. The presentation 
is based on a simple theoretical model (Neel 1949, Stoner 1948). In Section[3]wc 
give a brief overview of the most common magnetic characterization techniques 
and explain the information extracted from each one. In Section|3]we discuss the 
response of a dense MNP assembly to a magnetic field, when the interparticle 
interactions are important and lead to a collective behavior of the MNPs. Mean- 
field models are presented and an introduction to modern numerical techniques 
(Monte Carlo, Magnetization dynamics) to tackle this problem are presented. 
The chapter is summarized in Section [S] and the perspectives in this field are 
presented in Sections [51 

2 Isolated magnetic nanoparticles 

In this section we derive the criterion for formation of single-domain MNPs and 
examine the magnetization process at zero temperature by coherent rotation of 
magnetization (Stoner- Wohlfarth model). The behavior of a MNP assembly at 
finite temperature is discussed and the related concepts of supcrparamagnetism 
and blocking temperature are introduced. The effects of an applied dc magnetic 
field is examined within the simplest model assuming uniaxial anisotropy and 
bistability of particle moments (Neel model) . 

2.1 Single-Domain Particles 

The ground state magnetic structure of a ferromagnetic (FM) material is the 
outcome of the balance between three different types of energies, namely, the 
exchange (U ex ), the magnetostatic (U m ) and the anisotropy energy (U a ). The 
exchange interaction has its origin in the Pauli exclusion principle for electrons. 
Let the FM material be divided in small cubic elements each one carrying a 
magnetic moment jfc. The exchange interaction between the cubic elements 
favors parallel alignment of neighboring magnetic moments and it is written in 
the usual Heisenberg form as U ex = —(A/ a 2 ) J^ij cosOij, where A is the stiffness 
constant, a is the lattice constant and fty is the angle between moments at sites 
i and j. The stiffness constant is related to the microscopic exchange energy J 
through the relation A = zJS 2 /a, where S is the atomic spin and z = 1,2,4 
for sc, fee and bee lattice, respectively. The magnetostatic energy, is the sum of 
Coulomb energies between the magnetic moments comprising the FM material. 
It can be expressed as U m = —fioH c [ ■ M/2, where is the demagnetizing 
field and M the sample magnetization. The anisotropy energy, is the energy 
required to orient the magnetization at an angle (9) relative to certain fixed axes 
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of the system, known as the easy axes. The microscopic mechanisms leading 
to anisotropy can be quite diverse and the most common types of anisotropy 
found in FMs are as follows: 

(i) Crystal anisotropy. It arises from the combined effects of spin-orbit cou- 

pling and quenching of the orbital momentum that produce a preferred 
orientation of the magnetization along a symmetry axes of the underly- 
ing crystal. For a uniaxial materials (e.g. hexagonal Co) it has the form 
U a = K\siv?Q + K 2 sin 4 9 + where Ki,K 2 , ... are the anisotropy con- 
stants, and 9 the angle between the magnetization direction and the easy 
axis. Typical values for cobalt are Ki = 4.5 x 10 6 J/m 3 and K% = 
1.5 x 10 5 J/m 3 . For cubic crystals (e.g. fee Fe, Ni) it reads U a = 
K 1 (a 2 a 2 + a\a\ + a\a\) + K 2 a\a\a\ + where ai,a 2 ,as are the direc- 
tion cosines of the magnetization direction. Typical values for Fe are 
A'i = 4.8 x 10 4 J/m 3 and K 2 = ±0.5 x 10 4 J/m 3 . 

(ii) Stress anisotropy. It is produced by the presence of stress in the sample 
and it has a uniaxial character U a = K a sin 2 9, where K a = |Ai(T, with 
the magnetically induced isotropic strain and a the stress. 

(iii) Surface anisotropy. This is caused by the presence of sample free bound- 
aries, where the reduced symmetry and the presence of defects can induce 
additional anisotropy. It is important in MNPs because of the substantial 
surface-to- volume ratio. 

(iv) Shape anisotropy. This occurs because on one hand the demagnetizing 
field depends on the shape of the magnetized body and takes the low- 
est value along the longest axis of the sample, and on the other hand, 
U m is minimized when M is parallel to Hd- As an example, consider a 
specimen in the shape of prolate spheroid with major axis c and minor 
axis a, magnetized at an angle 6 with respect to c-axis. Then, U m = 
if [N C (M cos 6*) 2 + N a (M sm9) 2 ] = \{N C - N a )M 2 sin 2 9 , where N c and 
N a are the demagnetizing factors along the corresponding axes. This ex- 
pression for U m has the form of uniaxial anisotropy with K s = h(N c — 
N a )M 2 . Typical cases, are a spherical specimen with K s = 0, an infinitely 
thin planar specimen with iV|| = (in- plane) and N±_ = 1, and a infinitely 
long (needle-shaped) specimen with N\\ = ^ (along the axis) and N± = 0. 

In studies of the magnetic properties of MNPs, it is a common practice, to 
describe, within the simplest approximation, the overall effect of the various 
anisotropy types by an effective uniaxial anisotropy term U a — K e ff sin 2 9. 
The constant K e ff accounts for the total effect of crystalline, surface and shape 
anisotropy. 

A bulk FM material is composed of many uniformly magnetized regions 
(domains). The direction of magnetization in different domains varies, and in 
a bulk sample it is randomly distributed leading to a non-magnetized sample 
even at temperatures far below the Curie point. The formation of magnetic 
domains in FM materials results from the competition between the exchange and 
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Figure 1: Onc-dimcnsional model of a FM. (a) Long-range order, (b) An in- 
finitely thin DW (dashed line). The increase of exchange energy at the wall is 
higher than the decrease of the magnetostatic energy, (c) A 180° domain wall 
spread over N = 10 sites. The gradual rotation of atomic moments produces a 
state with lower total energy compared to (b). 



the magnetostatic energy. The former favors perfect alignment of neighboring 
moments and the latter is reduced by breaking a uniformly magnetized body 
into as many as possible regions with opposite magnetization directions. The 
outcome of this competition is the formation of a certain number of domains 
in a sample with a particular orientation of the magnetization directions. A 
typical domain size in a bulk ferromagnet is 1/im. 

Neighboring magnetic domains are separated by a region where the lo- 
cal magnetization changes gradually direction between the two opposite sides, 
known as domain wall (DW). Domain walls have finite width (5 W ) determined 
by the balance between the exchange and anisotropy energy. As an example, 
consider an one-dimensional model of a DW in a uniaxial material, where a 180° 
rotation of magnetization is distributed over N sites, as shown in Fig. [TJ The 
total energy per unit area reads 



For a typical exchange stiffness value (A 10~ n J/m), Eq.Q predicts for iron 
S w rj 0.4 fim while for a magnetically harder material like cobalt, 8 W ~ 60 nm. 
Substituting the result of Eq. @ into in Eq. ([T]) provides the areal energy density 
of the DW 



a(N) = a ex +o- a = JS 2 {n/N) 2 {N/a 2 ) + NaK^ 



(1) 



Minimization with respect to N leads to 



S w =Na = 7r(A//^i) 1/2 . 



(2) 



<7w = 



2v{AK 1 ) 1 ' 2 



(3) 
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Consider a finite sample of a FM material, with size d. As the size of the 
sample is reduced, the number of DWs it contains decreases, because fewer 
regions with opposite directions of magnetization are required to reduce the 
magnctostatic energy. Below a critical value of the system size, the sample does 
not contain any DW and it is in a single domain (SD) state exhibiting satu- 
ration magnetization (M s ). For a spherical particle, the critical diameter (d c ) 
can be estimated as follows: the SD state is stable when the energy needed to 
create a DW that spans the whole particle, U w = a w irr 2 , is greater than the 
magnetostatic energy gain from the reduction to a multidomain state, which is 
approximately equal to the magnetostatic energy stored in a uniformly mag- 
netized sphere, U m = ^fioM^V, with M s the saturation magnetization and 
V = 4^r 3 . The condition U w = U m provides 



(AK 



1 1/2 



(4) 



For Fe, this approximation gives r c sa 3 ran, which is by far too small. The rea- 
son is that the DW is assumed to have the same one-dimensional structure as in 
the bulk material. An improved calculation that considers a three-dimensional 
confinement of the DW provides for the critical radius: 



/ 9A 



I ^ ) 

a 



(5) 



In the case of Fe, numerical solution of Eq.([5]) gives r c w 25 ran, which is 
very close to more accurate micromagnetic calculations and the experimentally 
obtained value (Cullity 1972). 



2.2 Magnetization by Coherent Rotation 

The magnetization (M) of a bulk FM crystal that contains many magnetic do- 
mains, changes under application of an external magnetic field (H), a process 
known as technical magnetization. However, the value of M is not a unique 
function of H and the state of the sample prior to application of the field is 
important. This is the phenomenon of magnetic hysteresis, which is commonly 
depicted by drawing the M — H dependence under a cyclic variation of the field 
from a positive to a negative and back to a positive saturation value (hysteresis 
loop). Two important characteristic values of a hysteresis loop are the rema- 
nence (M r ), namely the magnetization after removal of the saturating field, and 
the coercivity (H c ), namely the field required for the magnetization to vanish. 
In a bulk FM crystal, the magnetization proceeds by two basic mechanisms, 
namely domain wall motion (weak fields) and rotation of magnetization (strong 
fields). 

In MNPs, the change of magnetization under an applied field proceeds only 
by rotation, because formation of DWs is energetically unfavorable. During the 
magnetization rotation the atomic moments of the MNP remain parallel to each 



other and the MNP behaves as a giant molecule carrying a magnetic moment 
of a few thousand Bohr magnetons (/x ~ 10 4 /is for a 5 nm diameter Fe MNP). 
This process of magnetization is known as coherent rotation or Stoner- Wohlfarth 
(SW) model, after the authors who introduced and solved it (Stoner 1948). We 
discuss it briefly next. Consider a MNP with uniaxial (effective) anisotropy K\ 
along an easy axis taken to be the z-axis (Fig. [5]). For an applied field that 
makes an angle (9o with the easy axis, we wish to determine the equilibrium 
position of the magnetic moment (i = M S V. Let ~jt make an angle 9 with the 
easy axis, then the total energy density reads 

u = -K x cos 2 (9 - 6> ) - H HM S cos 6 (6) 

The equilibrium condition (zero-torque) is 

du 

— = 2K lS in(6 - 9 ) cos(0 - 9 ) + /j, HM s sin6 = (7) 
dO 

and introducing the dimensionlcss quantity h = H/H a with the anisotropy field 
H a = 2Ki/fj, M s , Eq.© becomes 

sin(2(0-0 o )) + 2/jsin0 = O. (8) 

We define the reduced magnetization along the field m = \i cos 6 /M S V = cosQ 
and the solution of Eq. © is written as 

2m(l - TO 2 ) 1/2 cos26> + sin26» (l - 2?ti 2 ) + 2h{l - m 2 ) 1 ' 2 = (9) 

The remanencc (h — 0) and coercivity (m = 0) are readily obtained from Eq.@ 
as 

m r = cos^o and h c = sin cos ■ (10) 

For non-zero field values, Eq.© is solved for as a function of m and the data 
are shown in Fig. [2] Consider the two extreme cases, namely for 6*o = 90° 
(hard-axis magnetization) and 9q = 0° (easy-axis magnetization) . In the former 
case, the magnetization shows zero coercivity and a linear field dependence. 
In the latter case, the magnetization remains constant until the reversing field 
becomes equal to the anisotropy field, and then an irreversible jump of the 
reduced magnetization from m = +1 to m = —1 is seen. These extreme cases 
demonstrate the distinct mechanism of switching by rotation that can occur in 
an assembly. More generally, at an arbitrary field angle, an irreversible jump 
of the magnetization occurs at the so called switching field (H s ) defined as the 
field value satisfying dm/dh — > oo. At H = H s the local minimum of the total 
energy, corresponding to the higher energy state (magnetization opposite to the 
applied field) disappears and the system jumps to the remaining minimum that 
corresponds to a magnetization direction along the field (see Fig. [3]) . In other 
words, H s is an instability point of the total energy and it satisfies du/dQ = 
and d 2 u/d6 2 = 0. In the SW model, the stability condition reads 

d 2 u 

— r =O=>cos2(0-0 o )±/isin0 = O (11) 
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Figure 2: (a) Sketch of a magnetic nanoparticle with uniaxial anisotropy along 
the z-axis and an applied field at an angle (9q) with respect to the easy axis, 
(b) Magnetization curves within the Stoner-Wohlfarth model for various field 
directions. The initial direction of the magnetization is taken along the field. 



From Eqs.© and (TTT]) we obtain for the switching field h s = H s /H a 

h s = (cos 2 ^e + sm 2 ^e Q r 3 / 2 (12) 

By comparison of Eqs.(flO]l and ([12]) one finds that h c < h s for 45° < 9 < 90° , 
namely switching happens after the magnetization changes sign, while for field 
angles close to the easy axis, 0° < 6q < 45°, the magnetization changes sign 
by an irreversible jump (h c = h s ). The physical distinction between h s and 
h c can be understood by the following example. Consider a SW particle under 
application of a reversing field h = /i c , which brings the particle's moment ~jl in 
a direction perpendicular to the field, so that m = 0. Then the field is switched 
off adiabatically. If h c < h s (i.e. 45° < 9 < 90°), ~jt will return back to the 
positive rcmanence value (m = +1), while if h c — h s (i.e. 0° < 6>o < 45°), 
~ji will jump to the negative remanence state (m = —1). The switching field 
of a hard (i.e. large anisotropy) magnetic material is a physical quantity with 
great technological interest in magnetic recording applications. In these, the 
information bit is stored in the direction of magnetization and the switching 
field is the field required to write or erase this information. 

Stoner and Wohlfarth (Stoner 1948) also studied an assembly of isolated 
MNPs with easy axes directions distributed uniformly on a sphere [random 
anisotropy model, RIM). The reported values for the remanence and coercivity 
are 

m r = 0.5 and h c = 0.48. (13) 
This result is particularly useful as random easy axis distribution is found 
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Figure 3: Dependence of total energy on the direction of the particle's moment 
(see Eq.©), for various strengths of the applied field (h = H/H a ). The energy 
minimum at 6 = 7r becomes unstable at the switching field h s = 1. 



in most MNP-based materials (granular films, cluster-assembled films, self- 
assembled arrays, etc) 

As a final remark, we remind that in the SW model thermal effects are 
ignored (T = 0), thus energy- minimization with respect to the magnetic moment 
direction is a sufficient condition to determine the field-dependent magnetization 
at equilibrium. The magnetic behavior of SD particles at finite temperature is 
discussed in the following section. 

2.3 Magnetic behavior at finite temperature 

How do thermal fluctuations affect the average magnetization direction of an 
isolated MNP ? How does the presence of an applied field modify the magnetic 
response at finite temperature ? Is the assembly magnetization stable in time, 
when the MNP moment arc subject to thermal fluctuations ? These points are 
briefly discussed next, along the lines of a model first studied by Neel (Neel 
1949). 

2.3.1 Superparamagnetism and Blocking temperature 

Consider an assembly of identical SD particles with uniaxial anisotropy. The en- 
ergy (per particle) is U = — K{V cos 2 8, where 9 is the angle between the single 
particle magnetic moment ft and the easy axis. The energy barrier that must be 
overcome for a MNP to rotate its magnetization is Ef, = K{V. As first pointed 
out by Neel (Neel 1949), thermal fluctuations could provide the required energy 
to overcome the anisotropy barrier and spontaneously (i.e. without externally 
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applied field) reverse the magnetization of a MNP from one easy direction to the 
other. This phenomenon can be thought of as a Brownian motion of a particle's 
magnetic moment. The assembly shows paramagnetic behavior, however it is 
the giant moments of the MNPs that fluctuate rather than the atomic moments 
of a classical bulk paramagnetic material. This magnetic behavior of the MNPs 
is called superparamagnetism (SPM) (Bean 1959) At high enough temperature. 
ksT » K{V, the anisotropy energy can be neglected and the assembly magne- 
tization can be described by the well known Langevin function M = nM s L(x), 
where n is the particle number density, and x = i^lq^lH /ksT . Thus, the features 
serving as signature of superparamagnetism are the scaling of magnetization 
curves with H/T, as dictated by the Langevin function, and the lack of hys- 
teresis, i.e. vanishing remanence and coercivity. Moreover, the major difference 
between classical paramagnetism of bulk materials and SPM is the weak fields 
(H ~ 0.1 T) required to achieve saturation of a MNP assembly magnetization 
M. This occurs because of the large particle moment (fi ~ 10 4 /i_e) compared 
to the atomic moments (^ at ~ /Ub)- 

Measurement of magnetization curves at sufficiently high temperature can, 
in principle, be used to extract the particle moment /i. In practice, two compli- 
cation arise. First, the presence of different particle sizes in any sample produces 
a convolution of the Langevin function with the volume distribution function. 
Second, interparticle interactions, modify the reversal mechanism and the SW 
model needs extensions, which are discussed in the Section 

At low temperature, fcgT << K{V , the anisotropy barriers are very rarely 
overcome (weak thermal fluctuations), the assembly shows hysteresis and this 
is called the blocked state. 

One might now ask, whether there exists a temperature value that draws the 
border between the blocked and the SPM state. Following Neel's arguments, 
we assume that thermal activation over the anisotropy barrier can be described 
within the relaxation time approximation (or Arrhenius law) as 

t = T exp(KiV/k B T), (14) 

where l/2r is the probability per unit time for a reversal of ft . The intrinsic 
time To depends on the material parameters (magnetostriction constant, Young 
modulus, anisotropy constant and saturation magnetization). Typical values 
are To ~ 10 -10 — 10 _9 s as obtained by Neel. To detect the superparamagnetic 
behavior experimentally, the MNP must be probed for a long enough period of 
time to perform many switching events that would produce a vanishing small 
time-average magnetic moment. If T m is the measuring time-window, the con- 
dition for SPM behavior is r m 3> t. The strong (exponential) dependence of r 
on temperature (see Eq. (fl"4]) ) permits us to define a temperature value (or more 
precisely, a very narrow temperature range) above which the relaxation time is 
so small that SPM behavior is observed. This is called the blocking temperature 
(Tf,) of the assembly, and is given by 

T b = KiV/kBHTm/r). (15) 
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For T < Tb, the particle moments fluctuate without switching direction (on av- 
erage) and the assembly is in the blocked state exhibiting hysteresis. For T > T& 
the assembly is in the SPM state, hysteresis disappears and thermal equilibrium 
is established. It is remarkable, that the value of Tb depends on r m , which is a 
characteristic of the experimental technique adopted. For example, in dc sus- 
ceptibility measurements r m ps 100s, in ac susceptibility r m ps 10~ 8 — 10 4 s, 
in Mossbauer spectroscopy r m ps 10~ 9 — 10~ 7 s and in neutron spectroscopy 
r m ps 10~ 12 — 10~ 8 s. Therefore, if Tb is of interest for a particular application, 
the measurement technique implemented must imitate the real conditions. For 
example, to study the reliability of magnetic storage media, dc magnetic mea- 
surements over a wide time window (r m ~ 10 2 — 10 4 s) should be used, while 
to study magnetic recording speed, ac measurements are appropriate. 

Brown (Brown 1963) extended the treatment of thermal activation over the 
anisotropy barrier, allowing also for fluctuations of \i transverse to the easy axis, 
which Neel has neglected, and obtained a different expression for tq. However, 
the common feature of both studies is the temperature and volume dependence 
of t, so the final result, Eq. (Tn)) , is referred to as the Neel-Brown model. 

In a polydisperse assembly, the distribution of particle volumes /(V), pro- 
duces a corresponding distribution of blocking temperatures f(Tb). Then, at 
a certain temperature T the assembly contains a mixture of blocked and SPM 
particles. The MNPs with volumes above a critical value V c , fulfill the require- 
ment of strong thermal energy with respect to their anisotropy barrier, and are 
SPM, while those while those with V < V c are blocked. From Eq.([T4|), the 
critical volume reads V c = kbT \xi(r m / To) / 'K\. As explained above for Tb, also 
for V c the experimental determination depends on the technique adopted. Most 
preparation techniques result in polydisperse samples and the problem of ex- 
tracting the size distribution function from magnetic measurements, pioneered 
by Bean and Jacobs more than fifty years ago (Bean 1956) remains a difficult 
task mainly due to the complications introduced by interparticle interactions. 
Knobel and colleagues have recently reviewed this subject (Knobel 2008). 

2.3.2 Thermal relaxation under an applied field 

Consider an assembly of N identical MNPs with uniaxial anisotropy along 
the z-axis and let their moments point initially along the + z-axis. Assume 
that a magnetic field H, weaker than the switching field , which is equal to 
H a , is applied along the —z-axis. Then, the total energy per particle reads, 
U = — K\V cos 2 (# — 0q) + hqHM s V cosQ. It exhibits two non-equivalent local 
minima at 9 = 0,ir with values U± = —K\V ± M S VH and a maximum at 
9 = tt/2 with U m ax = KiV{H/H a ) 2 , as shown in Fig. H The energy barriers 
and the corresponding relaxation times for the forward (+) and the backward 
(— ) rotations are 

E b ± (H)=K 1 V(lTH/H a ) 2 and r± = r exp(E±/k B T). (16) 

The change of tq due to the field is much weaker than the change of the expo- 
nential factor and as such it is neglected in the above equation. 
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Figure 4: Total energy of an isolated particle with uniaxial anisotropy subject to 
a negative field parallel to the easy axis with value less than the switching field 
(0 < H < H s ). Energy barriers (Ef,) and relaxation times (r) for the forward 
(+) and the backward (— ) process are not equal. 



The blocking temperature, as measured within a time- window r m , is reached 
when the observation time equals the forward relaxation time r+, because the 
latter corresponds to a moment flip from the initial state along +z to the op- 
posite direction, namely a process that reduces the initial magnetization. From 
Eq. fTB)) one obtains 

T b (H) = Kl l { ] ~ H/H f ee T h (0)(l - HlH a f. (17) 

which indicates that the blocking temperature is reduced by the presence of 
a reverse field. By completely symmetric arguments one could show that Tj, 
increases in the presence of a field with the same direction as the initial magne- 
tization. 

Since thermal fluctuations act in synergy to a reverse field in switching the 
moment of a MNP, it is expected that the coercivity of an assembly will decay 
with temperature. As discussed above, for a particle with its moment along 
the +z-axis, a reverse field (0 < H < H a ) reduces the barrier for reversal to 
the value E£ given in Eq. ([T6|) . If the field is strong enough, it will reduce 
the barrier to the value appropriate for superparamagnetic relaxation, namely 
fc_BTln(r m /To), and the (time-average) magnetization will vanish. On the other 
hand, the reverse field that makes the magnetization vanish is by definition the 
coercive field. Therefore, the following relation holds 

K X V{1 - Ho/Ha) 2 = k B Tln{T m /T ) (18) 
which, using Eq. (|15[) . provides the temperature dependent coercivity 

H C (T) =H a \l- {T/n) 1 ' 2 ] . (19) 
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The microscopic mechanism of thermal activation of the MNP moment 
over the anisotropy barrier, produces a macroscopically measured time-decay 
of the magnetization. We derive this dependence assuming that when a mo- 
ment switches direction it continues to remain along the easy axis (Neel 1949). 
Then, at time t, N + particles occupy the lower minimum at 6 = 0, and the 
rest iV_ = N — N + particles occupy the higher minimum at 9 = it. The time- 
evolution of is governed by the rate equation 



dN + N + iV_ 

dt T + T_ 



(20) 



The magnetization per particle is given as M(t) = (2N + (t) / N — l)M s , and 
solution of Eq. (f20"|) provides 



M{t) = Moo + (M - Moo) exp(-t/r) (21) 
with 1/r = l/r+ + 1/t_ being the reduced relaxation time and 

Moo = T+ ~ T M s and M = ( ™ + }^ - 1 ) M s (22) 



t++t_ 3 u \N-(0) 

the time-asymptote and initial values of the particle magnetization, respectively. 
Eq. (|2ip indicates that the magnetization decays exponentially towards the equi- 
librium value Mqo, reached as t — > oo. In other words, equilibrium is reached 
when the population of the energy minima is proportional to the corresponding 
relaxation times (N+/N- = t+/t_), as dictated by Eq. ([20|) . When the applied 
field is strong enough (H > H s ) to produce only one minimum, thermal equilib- 
rium is always reached. Obviously, in the absence of an external field, thermal 
equilibrium is reached when the two equivalent minima are equally populated 
(iV_|_ = -/V_), resulting in a vanishing magnetization. 

Notice that in Eq. (f2"0f we assumed bistability of the moment direction, which 
is a valid approximation provided the anisotropy barrier is high (K±V 3> kgT). 
For lower anisotropy barriers or elevated temperature (K\V w fc^T), the trans- 
verse fluctuations of ~fi, or, in other words, intra- valley motion around the 
energy minimum should be taken into account. A general treatment of thermal 
relaxation of SD MNPs was pioneered by Brown (Brown 1963) and extended 
to the case of an applied external field (Aharoni 1965, Coffey 1998, Garannin 
1999). 

If an assembly is polydisperse, characterized by a volume distribution f(V), 
a distribution of blocking temperatures /(Tfc) exists. However, it remains still 
unclear if the mean value < Tt > is the appropriate blocking temperature of 
the assembly, which should be substituted, for example, in Eq. (|19p . This point 
is discussed further in the literature (Nunes et al 2004). 

In a polydisperse assembly, a distribution of relaxation times /(r) exists, 
with /(r)d(ln r) the probability of a MNP to have In r in the range (In t, In r + d In r) 
and the normalization condition J /(r)d(lnr) = 1. In this case the magneti- 
zation can be obtained by a superposition of the single-particle magnetization 
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properly weighted, as follows 

MO) = M s I [I- exp(-t/r)] ^dr, (23) 
Jo T 

where the term in brackets is the probability per unit time for a particle not to 
flip its moment. For a broad enough distribution, the observation time t will 
satisfy n <C t <C T2 , where t\ and T2 are the minimum and maximum relaxation 
times of the assembly, respectively. Assuming a uniform distribution /(r), it 
can be shown that the magnetization exhibits a logarithmic relaxation 

M(H, t) = M(H, 0) - S(H, T) ln(t/r ) (24) 

with S the magnetic viscosity of the system. Thus, polydispersity produces a 
much slower decay of magnetization with time. 

The discussion so far, refers to a field applied parallel to the easy axis. 
However, random anisotropy is most commonly found in MNP assemblies and 
the the necessity to study the effect of a tilted field with respect to the easy axis, 
arises. In this case, the calculation of the energy barriers and relaxation time 
is a much more complicated task and no analytical solution exists. Numerical 
studies (Pfeiffer 1990) showed that the energy barrier for an applied field at an 
angle 9q to the easy axis can be approximately written as 

E b (6 ) = K\V(1 - H /H a f m+1 - Xih > (25) 

where h s is given by Eq. (fT2")) . In the limit of 6 = 0, Eq. ([23| reduces to Eq. fTTJ)) . 

The temperature dependence of the coercivity for a monodisperse assembly 
with random anisotropy has also been obtained numerically (Pfeiffer 1990) as 

H C (T) = 0A8H a [1 - (T/Ti,) ' 77 ] , (26) 

which at T = reduces to the SW result of Eq. (fT3| . A detailed theoretical 
study of the relaxation time for a non-uniaxial applied field can be found in the 
review by Coffey and colleagues (Coffey 1993). 

As a concluding remark, the presence of polydispersity and random anisotropy 
makes the description of the magnetic behavior of an assembly intractable to 
exact analytical treatment. Instead, numerical approximations and simulation 
methods provide the alternative theoretical tools to study these systems. Nu- 
merical simulation approaches are introduced in Section 2J 

3 Magnetic Measurements 

Thermal relaxation has a dynamic character, therefore, the relation between the 
various relaxation times of the assembly and the measurement time is a decisive 
parameter for the outcome of a measurement. Additionally, if the assembly 
is not at equilibrium during the measurement, or if it changes its equilibrium 
state (for example, by adiabatic changes of the applied field) the result of the 
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measurements depends on the measurement protocol followed. In what follows 
we discuss two very common types of static measurements, that reveal the 
temperature and field dependence of the magnetization and provide evidence 
for superparamagnetic relaxation. Dynamic measurements are not discussed 
in this article. The interested reader can find more on the physical principles 
behind the most common magnetic measurement techniques in the review of 
Dormann and colleagues (Dormann 1997). 

3.1 Field-cooled (FC) and Zero-Field-Cooled (ZFC) Mag- 
netization 

This is a measurement protocol adopted for investigation of the temperature 
dependent magnetization of an assembly and it reveals superparamagnetic be- 
havior. It is performed in three stages. In the first, the sample is initially at 
a high enough temperature (T max ) to ensure a SPM state and it is cooled to 
low temperature (T m j n ) to approach its ground state. In the second stage, a 
weak field is applied (H <C H sat ), the sample is heated up to T max and the 
magnetization is measured as a function of temperature. This is the ZFC curve. 
In the third stage, the system is cooled down to T m j„, without removing the 
field, while the magnetization is recorded again, producing the FC curve. Dur- 
ing cooling and heating the temperature changes at the same constant rate. A 
typical ZFC-FC curve is shown in Fig. [5] As the temperature rises, the blocked 
magnetic moments align easier along the applied field leading to an initial in- 
crease of the ZFC curve. However, as soon as thermal fluctuations push the 
moments over the anisotropy barrier, thermal randomization of the moments 
produces a drop of the curve. Therefore, the peak of ZFC curve corresponds to 
the blocking temperature of the assembly. Notice that above T b the ZFC and 
FC curves coincide, because the system is in thermal equilibrium and the the 
heating (cooling) process is reversible. On cooling below T b the moments re- 
main partially aligned along the field, and the magnetization tends to a non-zero 
value. The magnetization vanishes at the ground state (T m i n ) if the measuring 
field is very weak, a random distribution of the easy axes exists and the as- 
sembly is non- interacting (dilute) . Deviations from any of the above conditions 
produce a non-zero value for M ZFC {T = 0). For isolated MNPs, the ZFC-FC 
curves are only weakly sensitive to the value of the applied field, provided that 
it is weak (H <C H s ). 

3.2 Remanent magnetization and Coercive field 

Remanent magnetization at a certain field, M r (H), is measured after switching 
off the previously applied field H . In an assembly of MNPs, remanence arises 
because the moments of some particles, which have rotated under an applied 
field and to do so they have overcome an energy barrier, cannot rotate back to 
their original direction after removal of the field. In a polydispcrse assembly, 
at finite temperature T, only the blocked MNPs, namely those with T& < T 
contribute to the remanence. Therefore, M r /M s = f™ f{E b )dE b where E b . c = 
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Figure 5: (a) Typical FC-ZFC magnetization curves. The curves join at the 
peak of the ZFC which corresponds to T&. Arrows indicate the direction the 
measurements are taken. For T > Tf, the system is in thermal equilibrium 
and the heating process is reversible, (b) FC-ZFC curves for a dilute (non- 
interacting) assembly of Fe nanoparticlcs (D — 3.0 nm,M s = 1720 emu/cc 
and K\ = 2.4 x 10 5 erg/cc). The blocking temperature (dotted line) decreases 
weakly with increasing measuring field. The non-zero values of Mzfc{T = 0) 
are due to the finite value of the measuring field. Data produced by Monte 
Carlo simulations (Section 14.31) 



K{V C is the critical barrier for SPM relaxation at temperature T. Taking into 
account that I& ~ V (see Eq. ([T5j) ). we deduce that 

dM r (T)/dT = f(T b ) (27) 

namely, the slope of M r (T) provides the barrier (or blocking temperature) dis- 
tribution function of the assembly. There are three different measurement proto- 
cols for the remanent magnetization, as first suggested by Wohlfarth (Wohlfarth 
1958) : 

(i) Thermoremanence TRM(H,T), measured at the end of a FC process with 

field H from T max down to the measuring temperature T. 

(ii) Isothermal Remanence IRM(H,T), measured at the end of ZFC process 
from T max down to the measuring temperature T, at which a field H is 
applied and then removed. 

(iii) DC Demagnetization remanence DcD(H,T). First, a ZFC process from 
T max down to the measuring temperature T is performed. Second, the 
sample is brought to saturation remanence IRM(oo,T). Third, a reverse 
field H is applied and then removed to leave the sample at the DcD(H, T) 
remanence. 
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Wohlfarth pointed out that for isolated MNPs the different remanent magne- 
tizations are related as DcD(H) = IRM(po) — 2 • IRM(H). More interestingly, 
the deviations from this equality, defined as 

AM (if) = DcD(H) - [IRM{oo) - 2 • IRM(H)} (28) 

quantify the character and strength of intcrparticlc interactions and arc obtained 
experimentally (O'Grady et al 1993). Positive AM values imply interactions 
with magnetizing character, and negative values indicate demagnetizing inter- 
actions. We should say that this is only a phenomcnological characterization 
of the interactions, because Eq.(|2"5f does not provide any information about 
their microscopic origin. However, Eq. ([2"5)l has been proved a standard tool 
for quantification of intcrparticlc interactions in complex MNP assemblies such 
as those used in modern industry of magnetic recording media (granular films, 
particulate media). Intcrparticle interactions are discussed in the Section [4] 

4 Interacting nanoparticle assemblies 
4.1 Introduction 

The magnetic interactions that are present in bulk magnetic materials pertain to 
MNP assemblies and they preserve their physical origin and characteristics. In 
particular, (direct) exchange between atomic moments separated by a few lat- 
tice constants can couple ferromagnetically or antiferromagnetically two MNPs 
via their surface atoms. Indirect exchange or Rudcrman-Kittel-Kasuya-Yosida 
(RKKY) interaction exists between MNPs hosted in a metallic matrix, which 
provides free electrons required to mediate the interaction between the atomic 
moments of the MNPs. Finally, magnetostatic interactions, which are of minor 
importance in bulk magnets due to their weakness, become the dominant in- 
teractions in MNP assemblies with well separated MNPs. This situation occurs 
for two reasons. First, the exchange interactions have a very short range (up to 
~ 5 A), the RKKY interactions have an oscillating FM/AFM character with a 
period of a few A, which renders to zero their average effect on the MNP vol- 
ume, so both have a weak effect in interparticle coupling. On the other hand, 
magnetostatic interactions, in the lowest approximation, namely the dipolar 
contribution, are proportional to the magnitude of the coupled magnetic mo- 
ments, which for SD particles has an enormously large value compared to the 
atomic moments ([Imnp ~ IO'Vb ~ 10 3 /i a tom)- 

Further on we discuss the effects of magnetostatic (dipolar) interactions on 
the magnetic properties of MNP assemblies and their interplay with single- 
particle anisotropy. The complexity of this problem arises from the long-range 
(~ 1/d 3 , d being the interparticle distance) and anisotropic character of the 
dipolar interactions, namely the dependence of interaction energy on the orien- 
tation of the moments relative to the bond joining the particle centers (Fig. |6]). 
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Figure 6: Ground state configuration of magnetic nanoparticles with elliptic 
shape coupled by magnetostatic (dipolar) forces. The easy axis coincides with 
the long axis of the ellipse (shape anisotropy). Dipolar coupling induces anti- 
ferromagnetic ordering when the moments are forced (by anisotropy) to remain 
normal to the bond, as in (a). FM ordering (nose-to-tail) is favored when the 
easy axes are parallel to the bond, as in (b). In an assembly with random 
anisotropy, as in (c), the moments are aligned along the local easy axes. Dipo- 
lar interactions have a complex effect, leading to misalignment of the moments 
with respect to the local easy axis. 



Understanding and controlling the effects of dipole-dipole interactions (DDI) 
in MNP assemblies is of paramount importance to modern technology of mag- 
netic recording media for two opposite reasons. First, DDI couple the MNPs 
of an assembly. The ultimate goal in magnetic recording applications is to ad- 
dress each MNP individually and treat it as a magnetic bit. In this case, DDI 
have a parasitic role and one wishes to estimate and reduce their impact in the 
magnetic properties of an assembly. On the contrary, magnetic logic devices, 
have been proposed and built that exploit the magnetostatic coupling between 
ordered MNP arrays (linear or planar) to transfer a magnetic bit (usually a 
flipped moment) between two distant points in the array (Cowburn 2006). In 
this case, DDI are of central importance and the goal is to enhance and tailor 
their effects. 

Over the last two decades, many research groups have prepared and mea- 
sured MNP assemblies in various forms (granular films, ferrofluids, cluster as- 
sembled films, self-assembled nanoparticles, lithographic arrays of magnetic 
dots) and studied the intrinsic factors (host and particle material, particle size, 
particle density) and the extrinsic factors (temperature, field, measurement pro- 
tocol) that control the magnetic behavior. In many of these studies the presence 
of magnetostatic interactions has been confirmed. Among the above mentioned 
systems, the self-assembled MNPs prepared by a synthetic route offer the ad- 
vantage of containing well separated MNPs with a very narrow size distribution 
(oy ~ 5 — 10%), so they are ideal systems to study DDI effects. Experimen- 
tal observations on self-assembled MNPs that have been attributed to DDI 
include, reduction of the remanence at low temperature (Held 2001), increase 
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of the blocking temperature (Murray 2001), increase of the barrier distribution 
width (Woods 2001), deviations of the zero-field cooled magnetization curves 
from the Curie behavior (Puntes 2001), and difference between the in-plane 
and out-of- plane remanence (Russier 2000). Long-range ferromagnetic order 
in linear chains (Russier 2003), and hexagonal arrays (Puntes 2004, Yamamoto 
2008) of dipolar coupled single-domain magnetic nanoparticlcs has been demon- 
strated, supporting the existence of a dipolar superferromagnetic ground state, 
characterized by ferromagnetic long-range order of the particle moments. 

Investigations of the static and dynamic magnetic properties of dipolar in- 
teracting nanoparticle assemblies brought up fundamental issues related to the 
existence of a ground state which shares common features with spin glasses, 
such as slow relaxation, memory and ageing effects (Sasaki 2005). The latter 
are magnetic systems characterized by disorder and competing interactions that 
produce an energy landscape with many local minima, considered responsible for 
the occurrence of these effects. Dipolar intcrparticle interactions in dense and 
random nanoparticle assemblies are believed to cause a spin-glass-like behavior 
(Dormann 1997). 

Theoretical models have been developed in an effort to explain these obser- 
vations and related previous ones in assemblies with randomly located MNPs 
(granular films, cluster-assembled films). On a microscopic level, the presence of 
DDI between MNPs modifies the magnetization switching mechanism, which for 
an isolated MNP obeys the Neel-Arrhenius model. When anisotropic MNPs are 
dipolar coupled, the reversal mechanism is determined by the interplay between 
the single-particle anisotropy energy (E a ~ K\V) and the dipolar interaction en- 
ergy (Ed ~ fHftj/iii)- For weak interactions (Ed <C E a ), the moments reverse 
independently by thermal activation over energy barriers, which are however 
modified due to DDI. This limiting case is treated within a mean-field approxi- 
mation and is discussed in Section [4721 For strong interactions (Ed ^> E a ), the 
single-particle reversal is no longer valid. Reversal of one particle can excite the 
reversal of others, and the assembly behaves in a collective manner. Many-body 
energy barriers exist in the system, with values that depend on the configuration 
of all moments. Their evaluation becomes a formidable task and numerical sim- 
ulations offer in this case an indispensable tool. Numerical methods are briefly 
discussed in Section l4~3l 

For a detailed review on the magnetic properties of dipolar interacting MNP 
assemblies the reader is referred to the relevant literature (Dormann 1997, Far- 
rell et al 2005, Knobel et al 2008, Kechrakos and Trohidou 2008). The role of 
magnctostatic interactions in patterned magnetic media has been reviewed by 
Martin et al (Martin 2003) 

4.2 Mean Field models 

In an early attempt to include the effect of interparticle interactions in the ther- 
mal relaxation of MNPs , Shtrikman and Wohlfarth (Shtrikman 1981) assumed 
that the single-particle anisotropy barrier of a MNP is increased by the Zeeman 
energy due to the interaction field Hi nt produced by the moments of neighbor- 
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ing particles. In this model, the Neel relaxation time is obtained from Eq. (fTB]) 
with the applied field H replaced by the interaction field Hi nt . The mean-field 
approximations consists in replacing Hi nt by its thermal average value, which 
in Neel's model is 



Hi nt = H int tanh(^o ^H mt /k B T) n fiH? nt /k B T (29) 

the latter approximation being valid for weak interaction fields. Substitution of 
Eq.(p9l into Eq. (jl6]) gives 



r cxp 



I^V + ^HlJkBT 



Using the approximation 1 + x « 1/(1 — x) we write Eq.pop in the form 

KiV 



twto exp 



k B (T-T Q ) 



(30) 



(31) 



with To = /io/i 2 Hf nt / ksKiV . Eq. (f3"Tj) . also known as the Vogel-Fulcher law, 
indicates that the relaxation time of an assembly of interacting MNPs is the 
same as that of the isolated MNPs at a lower temperature. 

In the Shtrikman-Wohlfarth model, the temperature To, or equivalently the 
thermal average Hf nt , is not related to the microscopic parameters of the as- 
sembly, i.e. particle location, and is treated as a phenomenological parameter, 
fitted to experimental data. Dormann and colleagues (Dormann 1988, Dor- 
mann 1997) developed a statistical model for the average barrier in a dipolar 
interacting assembly which quantifies the interaction field and provides for the 
single-particle energy barrier 

E b = K X V + ?iiaiA'/ s 2 VL(ai^ 2 /Vk B T) (32) 

with ri\ the number of nearest neighbors of a particle, a\ = x v /y/2, x v the 
volume concentration of the particles and L(-) the Langevin function. Eq. ([3"2"j) 
indicates that the anisotropy barrier is increased due to DDI, thus the model 
of Dormann et al predicts an increase of the blocking temperature due to DDI. 
This model behavior has been observed in almost all types of MNP assemblies, 
with a few exceptions (Hansen and M0rup 1998). 

More recently, Allia et al (Allia 2001) used also a phenomenological approach 
to describe a superparamagnetic assembly with weak DDI. Namely, an assembly 
in a regime that the remanence and coercivity vanish, but the field-dependent 
magnetization varies with concentration of MNPs indicating the presence of 
DDI. The authors (Allia 2001) suggested that the dipolar field changes at a high 
rate and in random direction and therefore acts similar to the thermal field. The 
effect is accounted for by an apparent increase of the system temperature. The 
magnetization at temperature T is given by M = M S L [fiH/ks{T + T*)], with 
T* related to the average dipolar energy via ksT* = ni/j, 2 /d 3 and obtained 
by a fitting procedure. This model interpreted successfully the magnetization 
behavior of Co nanoparticles in Cu matrix and established the existence of the 
interacting superparamagnet regime (Allia 2001). 
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4.3 Numerical Techniques 

The mean-field models have the advantage of providing analytical expressions 
suitable for extracting system parameters from the experimental data by a fit- 
ting process. However, they are not applicable to strongly dipolar systems and 
they do not account for collective effects. Numerical techniques on the other 
hand, have the major advantage that they treat rigorously the local and tempo- 
ral statistical fluctuations of the macroscopic quantities characterizing the MNP 
assembly and provide an efficient interpolation scheme between the weak and 
the strong interaction regimes. We discuss briefly two most common numerical 
approaches, the Monte Carlo (MC) method and the Magnetization Dynamics 
(MD) method. 

4.3.1 The Monte Carlo method 

Different algorithms that mimic thermal fluctuations of the degrees of freedom of 
a physical system by means of (pseudo)random numbers go under the umbrella 
of Monte Carlo techniques. In the case of MNPs, two widely used algorithms 
are the Metropolis Monte Carlo (MMC) and the Kinetic Monte Carlo (KMC). 
The former is appropriate for a description of the equilibrium behavior of an 
assembly, while the latter also accounts, within a certain time scale, for the tran- 
sition to equilibrium. Both algorithms provide thermal averages of macroscopic 
quantities of interest in the canonical ensemble, i.e. at constant temperature. 
To do so a sampling of the phase space is performed, however the sampling 
procedures differ, as outlined below. 

The MMC algorithm samples the phase space, visiting preferentially states 
close to the equilibrium states (Importance Sampling). This is achieved when 
subsequently visited states form a Markov chain, meaning that the probability of 
visiting the next state depends only on the last visited one. To do so, one chooses 
the transition from state s to s' to occur with certainty, if it reduces the total 
energy (E s > < E s ) and with a finite probability p(s — > s') = cxp(— Es ' kT Es ), if it 
increases the total energy (E s i > E s ). Thus the system is allowed to climb-up 
energy barriers and slide-down toward energy minima until it reaches eventually 
the global minimum. 

In KMC the system jumps from a state s at a local minimum to a new state s' 
being also a local minimum by overcoming a barrier Eh- The jump is performed 
within a predefined time step At with probability p(At) = 1 — exp(— At/r) 
where t is the corresponding relaxation time with Arrhenius behavior, r = 
T cME b /kT). 

In both algorithms, interparticle interactions are included by replacing the 
applied field H with the total field Hi = H + Yljtj&i) Hint,ij, which includes 
the contribution from the interaction field H int i j. In contrast to mean-field 
theories, in MC and MD (see next section) techniques the interaction field is 
treated exactly, meaning that its value depends on the configuration of all the 
moments of the assembly and it changes at each time-step. 

An important distinction between the MC algorithms is that KMC simu- 
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lates the relaxation of the system in physical time, while time quantification of 
the MMC time step is possible only in the absence of interparticle interactions 
(Nowak 2000, Chubykalo 2003). However, a serious difficulty in KMC arises 
from the calculation of the local energy barrier required to obtain the transi- 
tion probability. In an interacting system the barrier depends on all degrees of 
freedom and its calculation is a formidable task (Chubykalo 2004, Jensen 2006), 
usually performed in an approximate manner (Pfeiffer 1990, Chantrcll 2001). 
Furthermore, the KMC assumes that the system evolves through thermally ac- 
tivated jumps over energy barriers, an approximation that becomes invalid at 
elevated temperatures (kT ~ E^), or when collective behavior becomes impor- 
tant, as, for example, in strongly interacting MNPs. Collective effects are better 
described within the MMC algorithm. 

For a detailed description and technical implementation of MC algorithms 
the interested reader could refer to the book by Landau and Binder (Landau 
2000) 

4.3.2 The Magnetization Dynamics method 

In this method, the equations of motion for the magnetic moments are integrated 
in time and time averages of the macroscopic magnetization are recorded. At 
zero temperature, the time-evolution of a magnetic moment fii under a total 
field Hi is described by the Landau-Lifshitz-Gilbert (LLG) equation 

^ = -A(t x li) - Bill, x (ti x Hi) (33) 

with A = 7/(1 + a 2 ). Bi = aj/(l + a 2 )[ii, 7 the gyromagnetic ratio, and a a 
dimcnsionlcss damping parameter. The first term on the r.h.s. of Eq. (|33p is the 
torque term leading to precession around the field axis and the second one is 
a phenomcnological damping torque that tends to align the processing moment 
with the field Hi. 

The dynamics at finite temperature are described by introduction of an 
additional field (-ff/,i) term in Eq. ([33|) with stochastic character. Hf is assumed 
to have zero time-average (white noise) and its values at different sites i,j or 
different instants t, t' are uncorrelated. The LLG equation augmented by the 
thermal field term is commonly referred to as the Langevin or stochastic LLG 
equation. 

4.3.3 Time scale of numerical methods 

The MC and MD techniques are complementary since they describe thermal 
relaxation of magnetic properties in different time scales. In the MD method, 
the characteristic time is a fraction (~ 10 -2 ) of the precessional (Larmor) period 
(~ 10~ 10 s), which implies that simulation times up to ~ 1 ns are presently 
attainable. Thus, MD is the appropriate scheme to investigate fast-relaxation 
phenomena as, for example, the reversal path of magnetization under an applied 
short field pulse (Berkov 2002, Suess 2002). 
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In KMC the characteristic time is the single-particle relaxation time (see 
Eq. (fT4|) '). which is much larger than To ~ 10~ 10 s in the temperature range 
of interest (kT <C Ef,), a fact that makes the method suitable to treat slow- 
relaxation problems, as, for example the thermal decay of magnetization in 
permanent magnets, a phenomenon that evolves within days or years (Van de 
Veerdonk 2002). 

Finally, when static magnetic properties are concerned, the system is at a 
stable (or metastable) state and the MMC is a powerful and sufficient scheme to 
describe, for example, long range order at the ground state or collective behavior 
at finite temperature (Kechrakos 1998, Jensen 2003). 

4.3.4 MMC study of dipolar interacting assemblies : A case study 

In this section we show typical results from MMC simulations of the magnetic 
properties of dipolar interacting MNP assemblies (Kechrakos 1998, Kechrakos 
2002). Our system contains N identical SD NPs with diameter D and uniaxial 
anisotropy in a random direction. The MNPs are located randomly in space or 
on the vertices of a hexagonal lattice. The former is an appropriate model for 
granular samples, and the latter for self-assembled MNPs. The total energy of 
the system is 

E = g J2 ^ % = :{ ' 5 / , ! /l '" il - S ' ' - k ]T(g, ■ el) 2 - h ]T(5 a ■ H) (34) 

ij l 3 i i 

where Si is the magnetic moment direction (spin) of the z-th particle, ei is the 
easy-axis direction, and Rij is the center-to-center distance between particles 
i and j. Hats indicate unit vectors. The energy parameters entering Eq.(jMJ) 
are: (i) the dipolar energy g = ^ofi 2 /Aird 3 , with /i = M S V the particle moment 
and d the minimum intcrparticlc distance, (ii) the anisotropy energy k = K\V , 
and (iii) the Zeeman energy h = Ho^H due to the applied dc field H. The 
energy parameters (g,k,h) entering Eq. (|34p . the thermal energy /csT, and the 
treatment history of the sample determine the micromagnetic configuration at a 
certain temperature and field. The freedom to choose an arbitrary energy scale 
makes the numerical results applicable to a class of materials with the same 
parameter ratios rather than to a specific material. The crucial parameter that 
determines the transition from single-particle to collective behavior is the ratio 
of the dipolar to the anisotropy energy (g/k). 

We show in Fig. [7] the concentration and temperature dependence of the 
rcmancncc magnetization of a random assembly. Notice in Fig. [7^,, that weak 
DDI produce an increase of the remanence with concentration, while strong DDI 
have the opposite effect. Remarkably, the presence of free sample boundaries, 
can reverse the increasing trend of the remanence, due to the presence of a de- 
magnetizing field. When DDI are much stronger than single-particle anisotropy 
(g/k ~ 10), the remanence value is sensitive to the morphology of the assem- 
bly, as the peak around the percolation threshold indicates. This behavior is 
explained by the anisotropic character of DDI (see Fig. [5]). In Fig. [TJd, DDI 
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Figure 7: Dependence of the saturation isothermal remanence of a random 
assembly (a) on volume fraction of MNPs, at very low temperature (t/k=0.001) 
and, (b) on temperature, for fixed volume fraction. The particles have random 
anisotropy. The data are obtained by MMC simulations. 



interactions are shown to produce a much slower temperature decay, producing 
finite remanence values above the blocking temperature of the isolated MNPs. 
This result supports the predictions of the mean-field theory about the increase 
of the measured blocking temperature in dipolar interacting systems (Dormann 
1988). 

In chemically-prepared, self-assembled MNPs the possibility to control the 
interparticle separation by variation of the surfactant (Willard 2004) offers the 
possibility to study the dependence of on interparticle spacing while pre- 
serving the geometrical arrangement of the assembly (hexagonal). In Fig. [5] 
we show results for the ZFC magnetization (Mzfc) an d the blocking tempera- 
ture as obtained from the peak of the Mzfc(T) curve for a hexagonal array of 
dipolar interacting MNPs with random anisotropy. Parameters corresponding 
to Co nanoparticlcs are used (Kcchrakos 2002). The characteristic dependence 
of T& on the inverse cube of interparticle spacing can be used as a proof of the 
dominant character of DDI in an assembly. Notice also that Mzfc(T « 0) 
assumes a positive value that increases with d values. This feature arises from 
the gradual formation of a long range ferromagnetic ground state, due to DDI. 

More examples of MC or MD simulations and comparison to experiments 
on MNP assemblies can be found in the relevant literature (Vedmedenko 2007, 
Kechrakos and Trohidou 2008). 
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Figure 8: ZFC curves and blocking temperature for an ordered (hexagonal) 
assembly of identical MNPs with diameter D and center-to-center distance d. 
(a) Evolution of ZFC curves with decreasing d values, (b) Linear scaling of T& 
with inverse cube of d. Data obtained by MMC simulations. 



5 Summary 

We have discussed the main theoretical concepts that pertain to the magnetiza- 
tion properties of isolated (non-interacting) nanoparticles and their assemblies. 
We estimated the critical radius for formation of single domain particles and 
studied the zero-temperature magnetization reversal mechanism of coherent ro- 
tation (Stoner-Wohlfarth model), the thermally activated reversal (Neel-Brown 
model) and the related phenomenon of supcrparamagnctism occurring above the 
blocking temperature. Complications arising from size polydispcrsity, distribu- 
tion of easy axes directions, and applied field on the relaxation time for mag- 
netization reversal were discussed. Two standard experimental techniques for 
(static) magnetic measurements, namely the field and temperature dependence 
of magnetization were outlined. Finally, the subject of interparticle dipolar in- 
teractions was introduced along with the most common theoretical techniques 
used to analyze interacting systems. Examples from Monte Carlo studies of 
MNP assemblies were given. 

6 Future perspectives 

The dynamic behavior of MNPs in the presence of interparticle interactions is 
expected to remain a topic of intense scientific and technological research in 
the coming years. The research effort is expected to focus on both the atomic 
scale properties of individual magnetic nanoparticles and on the mcsoscopic 
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properties of nanoparticle assemblies. 

On the atomic scale, future goals will include : (i) Reduction of the mag- 
netic particle size without violating thermal stability (superparamagnetic limit) 
at room temperature . The technological benefit from progress in this direction 
will be the development of magnetic data-storage media with higher areal den- 
sity. Given that the SPM effect is not observed below a certain size of a MNP. 
due to disorder effects on the particle surface, the search for new high-anisotropy 
materials is required. Composite nanoparticles with a core-shell morphology 
(Skumryev 2003) constitute an interesting perspective. 

(ii) Understanding and control of surface effects. With reduction of particle size 
the contribution from surface moments become of increasing importance. The 
chemical structure of the surface (disorder, defects) controls the magnitude and 
type of the surface anisotropy, which is usually much (up to ~ 10 times) larger 
than the core anisotropy. Synthetic methods can offer indispensable routes to 
surface structure modification. Ab-initio electronic structure calculations are a 
valuable tool to predict the surface anisotropy values and modeling of MNPs as 
multi-spin system will reveal complex magnetization reversal mechanisms be- 
yond the Stoner-Wohlfarth model (Kachkachi 2000). Experiments on individual 
nanoparticles (Wernsdorfer 2000) offer a unique test of the above theories. 

The future task on the mesoscopic scale will be to understand and control col- 
lective magnetic behavior in ordered nanostructures (self-assembled MNPs and 
magnetic patterned media). Ordered nanostructures include chemically pre- 
pared self-assembled MNPs and lithographically prepared magnetic patterned 
media. Chemical synthesis of MNPs and self-assembly (bottom-up approach) 
is a very promising and cost-effective method to produce ordered MNP arrays 
(WiUard 2004). However, deeper understanding and improvement of the self- 
assembly process is required in order to achieve larger (beyond 1mm 2 ) sample 
area with structural coherence. There is still a remaining problem as nanopar- 
ticles self-assemble into hexagonal arrays that are incompatible with the square 
arrangements required in industrial applications. A resolution to this problem 
could be the recently demonstrated templated assembly (Cheng 2004). Litho- 
graphic patterning (top-down approach) offers better control over the geomet- 
rical aspects of the assembly but cannot yet produce nanostructures with size 
below ~ 10071TO (Martin 2003). Increase of lithographic resolution is demanded 
in order to achieve patterned media with smaller (below 0.1/j.m) characteristic 
size. On the measurements side, improvement of existing techniques to probe 
mesoscopic magnetic order and excitations is demanded. Recent examples are 
the observation of mesoscopic sale magnetic order in self-assembled Co nanopar- 
ticles by an indirect method (small-angle neutron scattering) (Sachan 2008) and 
by direct methods such as magnetic force microscopy (Puntes 2004) and electron 
holography (Yamamoto 2008). From the point of view of basic physics, ordered 
nanostructures constitute model systems to study collective magnetic behavior 
driven by magnetostatic interactions, because the size, the shape and the spatial 
arrangement of the magnetic nanostructures is well controlled. Known phenom- 
ena are to be demonstrated on the mesoscopic scale and new ones possibly to 
be discovered. As a recent example, we refer to the observation of magnetic 
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frustration is magnetostatically coupled magnetic microrods (Wang 2006), a 
phenomenon previously met in bulk magnetic random alloys (spin glasses). 

Finally, progress in numerical modeling will provide methods for bridging the 
atomic scale and the mesoscopic scale simulations. Such multi-scale simulations 
point to the future of theoretical investigations in the field of relaxation in 
magnetic nanoparticles and have only recently started to appear (Yanes 2007, 
Kazantscva 2008). 
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